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ABSTRACT 

bJO. 

^ ' The flat spectrum radio quasar PKS 1222+216 (4C+21.35, z = 0.432) was detected in the very high energy y-ray band by MAGIC 

■^C , during a highly active y-ray phase following an alert by the Large Area Telescope (LAT) onboard Fermi. Its relatively hard spectrum 

(70-400 GeV photon index F = 2.7 + 0.3) without a cut off, together with the observed variability on timescale of ~ 10 min challenges 
standard emission models. In particular, if the emission originates in a portion of the relativistic jet located inside the broad line region 
(BLR), severe absorption of y rays above few tens of GeV is expected due to the 77 e* process. These observations therefore 
imply the existence of a very compact (/?(, ~ 5 x 10'"* cm) and very fast blob located far beyond the BLR radius (to avoid absorption), 
responsible for the rapidly varying high energy flux. However the long term (days-weeks) coherent evolution of the GeV flux recorded 
I I by LAT indicates that there could be also the substantial contribution from another, larger, emission region. We model the spectral 

^JL^ ' energy distribution of PKS 1222+216 during the epoch of the MAGIC detection assuming three different scenarios, namely: (1) a 

one-zone model considering only the emission from a compact blob outside the BLR; (2) a two-zone model considering the compact 
blob plus an emitting region encompassing the whole jet cross-section located outside the BLR and (3) a two zone model with the 
jet emitting region inside the BLR. We find that in all cases the high-energy emission from the compact blob is dominated by the 
^ ' inverse Compton scattering of the IR thermal radiation of the dusty torus. Furthermore, both regions are matter dominated, with the 

, Poynting flux providing a negligible contribution to the total jet power. These results disfavor models in which the compact blob 

. is the result of^ reconnection events inside the jet or "needles" of high-energy electrons accelerated close to the BH. Instead, the 

' observational framework and our radiative models could be compatible with scenarios in which the jet is re-collimated and focussed 

at large distances from the central BH. 

' Key words, radiation mechanisms: non-thermal — 7-rays: galaxies — quasars: individual: 4C+21.35 
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1. Introduction this case inside the BLR, i.e. at distances d < /?blr, the en- 



O 



vironment is rather rich of optical-UV photons, ideal targets 

O Blazars emitting y rays at Very High Energy (VHE, E > 100 ^j. JC scattering with the relativistic electrons in the jet. If 

^ GeV) are still a mmority of the whole population, but their num- high-energy hump of FSRQs is dominated by this mecha- 

O ber IS steadily increasing due to the new generation of sen- ^^^^ (External Compton scenarios, EC, e.g., Sikora et al. 1994, 

^ ■ sitive Cherenkov aiTays (e.g. De Angelis et al. 2008, Hinton Ghisellini et al. 1998, Dermer et al. 2009) a pronounced soft- 

T-H , & Hofmann 2009). The majority of them are BL Lac objects, ^^j^g spectrum is robustly expected above few tens of 

^ . in particular those with the synchrotron bump peaking in the Q^y ^^^g ^^e decreased scattering efficiency (e.g. Tavecchio 

kI) ' UV/X-ray band (Highly peaked BL Lac, HBLs, Padovani & & Ghisellini 2008, Ghisellini & Tavecchio 2009). More impor- 

p><; Giommi 1995). The characteristic lack of important thermal ^^^^^ ^^^^^ radiation field makes the environment strongly 

H components in their spectra suggests that the jet propagates in ^^^^^^ ^ rays above few tens of GeV (e.g. Liu & Bai 2006, 

>y . a relatively "clean" environment. In this case it is widely be- Tavecchio & Mazin 2009, Poutanen & Stern 2010, but see Stern 

heved that the high energy component is due to the inverse ^ Poutanen 201 1). The combination of the two mechanisms, ef- 

Compton (IC) scattering of the synchrotron photons themselves fg^^j^g ^^^^^ ^jj^^jj^ thresholds, makes the detection of VHE 

(synchrotron self-Compton model, SSC, e.g. Maraschi et al. emission by FSRQs rather surprising. Even if the BLR is as- 

1992, Bloom & Marscher 1996, Tavecchio et al. 1998). ^^^^^ be "flat", as suggested in the past (e.g. Shields 1978) 

Pai-ticulai-ly interesting is the detection in the VHE band and by recent studies (e.g. Jarvis & McLure 2006, Decarii et 

of few blazars belonging to the class of Flat Spectrum Radio ^1. 201 1), the absorption cannot be avoided (Tavecchio et al., in 

Quasar (FSRQs). The broad emission lines observed in their prep) 
optical-UV spectra, typical of quasars, are though to originate in 

the so called broad line region (BLR), where dense (n ^ 10^"'° The prototypical FSRQ 3C 279 has been detected by 
cm"^) clouds of gas orbiting at typical velocities of few thou- MAGIC in two occasions, few days after optical flares (Albert 
sands km/s at distances 7?blr ~ lO'^""^ cm from the central et al. 2008, Aleksic et al. 2011a), with an extremely soft spec- 
black hole (BH) (e.g. Kaspi et al. 2007, Bentz et al. 2009) are trum and a flux variable on a timescale of a day. PKS 1510- 
photoionized by the UV continuum of the accretion disk. In 089 is another FSRQ detected at VHE by H.E.S.S. (Wagner & 
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Behera 2010). The recent detection of the FSRQ PKS 1222+216 
(also known as 4C+21.35,z = 0.431,Osterbrock&Pogge 1987) 
by MAGIC (Aleksic et al. 2011b) is particularly challenging. 
The VHE spectrum smoothly connects to the LAT spectrum 
(Tanaka et al. 2011) and it is well described by a rather hard 
power law, difficult to reconcile with the expected effects of 
absorption and decreased scattering efficiency discussed above. 
Moreover, the rapidly varying VHE flux, with a doubling of 
the flux in about 10 min, constrains the emitting region size to 
R < cfvar(l + z)6 - 2.5 X lO'^^CtJ/lO) cm (6 is flie relativistic 
Doppler factor). If these dimensions are related to the jet cross 
sectional radius, as in the standard one-zone scenarios, this im- 
plies that the emission region is located close to black hole, 
hence deeply inside the opaque BLR radiation field. 

As already proposed in Aleksic et al. (20 1 lb), a possible way 
out to the problem is to assume that the region producing VHE 
y rays lies beyond the BLR, at distances d > Rblr- At these dis- 
tances the extemal radiation field is likely dominated by the ther- 
mal radiation of the dusty torus reprocessing the disk radiation 
(e.g. Blazejowski et al. 2000) and both the opacity and the re- 
duced scattering efficiency problems are relaxed. However, even 
for extremely small values of the jet collimation angle (< 1 deg), 
at these distances we expect that the jet cross sectional radius 
is far exceeding the limit imposed by the rapid variability (see 
also Tavecchio et al. 2010, Foschini et al. 201 1). Hence we have 
to suppose either (as proposed by Ghisellini & Tavecchio 2008, 
Giannios et al. 2009, Marscher & Jorstad 2010) the existence of 
very compact emission regions embedded in the flow (as already 
done to explain the ultra-fast variability events observed in the 
two BL Lacs PKS 2155-304, Aharonian et al. 2007, and Mkn 
501, Albert et al. 2007) or strong recollimation and focusing of 
the flow (e.g. Stawarz et al. 2006, Bromberg & Levinson 2009, 
Nalewajko & Sikora 2009). 

It is interesting to note here that the problems posed by the 
interpretation of these observational evidences from blazars have 
strict analogies with the issues raised by the unexpected recent 
discovery of fast (~day) variability of the y-ray emission of the 
Crab nebula by AGILE (Tavani et al. 201 1) and Fermi (Abdo et 
al. 2011). In fact, as in blazars, in order to account for these 
short timescales it seems necessary to admit the existence of 
very compact knots inside the large scale relativistic wind from 
the pulsar (Vittorini et al. 201 1, Bednarek & Idek 201 1) whose 
non-thermal emission is possibly boosted by relativistic eflfects 
(Komissarov & Lyutikov 2011). It is therefore tempting to as- 
sume that similar physical processes, possibly related to plasma 
instabihties in relativistic flows, operate and shape the emission 
in rather difl'erent y-ray sources. 

In this paper we consider the problems posed by the obser- 
vation of PKS 1222+216 focusing on the question if standard 
leptonic models can reproduce the observed spectral energy dis- 
tribution (SED) and variabihty and, in particular, if a one-zone 
model is still a viable solution. To this aim we consider if a sin- 
gle, extremely compact emitting region (a "blob") located be- 
yond the BLR can account for the spectral and variability ob- 
served properties, concluding that, on the basis of the results 
of the SED modeling and on the power request, this solution 
could be viable. Alternatively, we study the scenario based on 
the existence of a two emitting regions, a "standard" emission 
region encompassing the whole jet cross section plus a compact 
"blob" responsible for the rapidly varying high-energy emission 
detected by MAGIC. This possibility is in part motivated by the 
observed long-term smooth evolution shown by the LAT light- 
curve, which is difficult to explain if the emission originates from 
uncorrected, fastly-evolving blobs. Moreover, we will show that 



this case is also energetically less demanding than the single blob 
model. This two-zone scenario resembles that already discussed 
by GhiselUni & Tavecchio (2008) for the case of PKS 2155- 
304. However, the existence of an intense extemal radiation field 
makes the situation somewhat different, because the radiative in- 
terplay between the two regions is less important than the EC 
emission with the IR photons (see Appendix B). 

We construct the SED using nearly simultaneous Swift 
(UVOT and XRT), Fermi/LAT and MAGIC data (§2). Then we 
apply the models to reproduce the SED (§3). Finally we discuss 
flie results (§4). 

We adopt J/o = 70 km s'^ Mpc-\ £Ia = 0.7, Qm = 0.3. 

2. The spectral energy distribution 

MAGIC detected PKS 1222+216 with a short (~ 30 min) obser- 
vation on 2010 June 17 (MJD 55364.9), during a period of ex- 
ceptional GeV activity (Tanaka et al. 2011). Data corrected for 
absorption by the interaction with the extragalactic background 
light modeled as in Dominguez et al. (201 1) are reported in Fig. 
1 (taken from Aleksic et al. 2011b). Aleksic et al. (2011b) also 
reported LAT data-points and the corresponding spectral "bow- 
tie" obtained for a short time of 2.5 hours encompassing the 
MAGIC observation. For comparison, the black thick solid line 
(from Tanaka et al. 2011) shows the LAT spectrum during qui- 
escence. 

Swift observed PKS 1222+216 several times in 2010 May- 
June, in correspondence with high activity of the source. 
Although there are no pointings exactly in time with the MAGIC 
observation, there is an observation few days after, on June 20, 
when the daily averaged LAT flux was similar to that on June 17 
(see the LAT light curve in Tanaka et al. 201 1). 

XRT (total exposure time of 4.5 ksec) data obtained dur- 
ing this observation (in PC mode) were analysed using the 
xrtpipeline task considering grade 0-12 events and using 
v6.10 version of the HEASOFT package and the calibration files 
caldb 9/2/2011. As commonly observed in FSRQs, the XRT 
spectrum (red) is relatively hard with a photon index F = 1.60 ± 
0.13 (the fit gives = 9.25 for 11 d.o.f.). The 0.3-10 keV flux 
is Fx = 5.5 X 10"'^ erg cm"^ s"'. For comparison we also ana- 
lyzed the data corresponding to a previous observation on May 
29. The flux was a factor « 1.7 higher. Fx = 8.3 x 10 erg 
cm~^ s"\ and the spectrum (cyan) was softer than on June 20, 
wifli F = 2.12 + 0.14 (x^ = 10.6 for 11 d.o.f.). It is interesting 
to note that while X-ray spectra of powerful FSRQs are gener- 
aUy hard, F < 1.5 (e.g. Kubo et al. 1998, Tavecchio et al. 2000, 
2002), soft X-ray slopes (indicating a dominant contribution of 
the synchrotron or SSC component in this band) are occasionally 
shown also by the other two TeV detected FSRQs, 3C279 (e.g. 
Aleksic et al. 201 la) and 1510-089 (e.g. Tavecchio et al. 2000). 

UVOT data-points (obtained with the standard procedure, 
e.g. Bonnoli et al. 2011) reveal a rather hard optical-UV spec- 
trum suggestive of the presence of direct thermal emission from 
the accretion disc (see e.g. Ghisellini et al. 2009, 201 1 for a dis- 
cussion). To further investigate this point we also consider SDSS 
data (taken on January 2008). To directly compare the UVOT 
and SDSS photometry we convert SDSS ugriz magnitudes into 
the standard Johnson system using the formulae provided by 
Chonis & Gaskell (2008). Using standard zero-points we fi- 
nally convert magnitudes into fluxes (magenta open squares). 
The SDSS data confirm the presence of a rather hard contin- 
uum. However, for the overlapping filters (U, B and V) there is 
a systematic diff'erence of about 0.8 magnitudes between UVOT 
and SDSS, translating into UVOT fluxes larger by a factor of 2, 
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Fig. 1. Spectral energy distribution of PKS 1222+216 close to 
the epoch of the MAGIC detection (2010 June 17). Red points at 
optical-UV and X-ray frequencies are from a Swift observation 
of June 20. For comparison, cyan data-points show the X-ray 
spectrum two weeks before, on May 29 (see text). Fermi/LAI 
(red squares and "bow tie") and MAGIC data (corrected for ab- 
sorption by the EBL using the model of Dominguez et al. 201 1) 
are taken from Aleksic et al. (2011b). The thick black solid 
line shows the LAT spectrum in quiescence (from Tanaka et 
al. 201 1). Magenta open squares are SDSS photometric points. 
Magenta filled pentagons are IR data from Malmrose et al. 
(2011). Green points report historical data (from NED, circles, 
and Tornikoski et al. 1996, squares). 



as clearly visible in the SEDs. This difference could reveal an in- 
crease of the accretion luminosity between the SDSS (Jan. 2008) 
and the UVOT (June 2010) observations, possibly related to the 
high activity in y rays. 

Malmrose et al. (2011) recently reported Spitzer observa- 
tions in the IR band for four sources, including PKS 1222+216. 
The IR data points (filled magenta pentagons in FigO) track a 
bump around 3 jjm which is well fitted by a black body with 
temperature of T = 1200 K, clearly related to the theimal emis- 
sion from the putative dusty torus. Finally, we also add historical 
radio (from NED) and millimeter (Tornikoski et al. 1996) data 
(green open circles and open squares, respectively). 

3. Modelling the SED 

3.1. Observational facts and problems 

In modeling the observed SED we are constrained/guided by the 
following observational facts: 

1) The MAGIC VHE spectrum (70-400 GeV) is well described 
by a hard power law, with photon index (after correction for 
absorption by the interaction with the extragalactic background 
light) of 2.7 + 0.3, and a cut-off for energies lower than 130 GeV 
is excluded. This spectrum smoothly connects with the LAT 
spectrum close to the MAGIC detection (Tanaka et al. 2011), 
strongly suggesting that high-energy and VHE emissions belong 
to a unique spectral component, originating in the same region. 

2) The MAGIC lightcurve shows a significative increase of the 
flux during the 30 min observation, with a doubling time of 



about fvar - 10 minutes. The causality relation R < cfvar(l + z)6 
allows us to constrain the size of the emitting region to R < 
2.5 X 10'"^(5/10) cm for typical values of the Doppler factor 
6^ 10. 

3) The LAT long-term lightcurve (Tanaka et al. 201 1) is charac- 
terized by periods of quiescence and smooth, long lasting (~ 1 
week) flares. The MAGIC detection coincides with the raising 
part of a flare lasting for approximately 3 days. The y-ray LAT 
flux (F>iooMeV ~ 6.5 X 10"^ ph cm"^ s"') was about half that 
recorded at the maximum of the flare (F>iooMeV ~ 13.5 x 10"^ 



>100MeV 

ph cm s '), reached the day after the MAGIC detection . 

Standard one-zone models for FSRQ generally assume that 
a single region in the jet, with a size comparable with that of the 
jet cross sectional radius, is responsible for the emission from 
IR to GeV frequencies. The location of this region is generally 
assumed to be inside the BLR (e.g. Dermer et al. 2009, Ghisellini 
& Tavecchio 2009), but scenarios considering regions beyond it 
have been discussed (e.g., Sikora et al. 2008, Marscher et al. 
2008). 

The observational facts listed above already pose some prob- 
lems to this view. Points 1) and 2) imply that the entire MeV- 
GeV and VHE emission component at the epoch of the MAGIC 
detection was produced in a very compact emission region out- 
side the BLR, to minimize the expected severe absorption above 
10 GeV (but see Stern & Poutanen 2011). In the framework of 
one-zone models, a first possibility is therefore to assume that 
the entire y-ray activity is due to the cumulative emission of very 
compact, uncorrelated traveling regions (resulting from, e.g. in- 
ternal shocks, Spada et al. 2001). However, in this case the ex- 
pected erratic behavior is in contrast with the smooth long-term 
evolution shown by LAT. One way to reconcile this scenario 
with point 3) seems to assume the existence of a very com- 
pact and stationary region: this would allow fast variations of 
the flux and, at the same time, the long term modulation of the 
jet power would account for the smooth and coherent evolution. 
As an alternative we could envision the existence of two emit- 
ting regions, a large region responsible for the long-term evolu- 
tion visible in the LAT band and an extremely compact region 
accounting for the fast variations. 

Motivated by the arguments above, in the following we 
present three different scenarios for the VHE flare of PKS 
1222+216 (see Fig.O. In the first case (A) we assume that the 
entire SED is produced by a single compact blob outside the 
BLR. In the other two cases we consider a two-zone model with 
the large region located outside (B) or inside (C) the BLR. For 
consistency with the scenario sketched above, in cases B and C 
we admit that the large region could substantially contribute to 
(even if not dominate) the LAT band also at the epoch of the 
MAGIC detection. 



3.2. Model setup 

A sketch of the assumed geometry is shown in Fig.|2] In all cases 
a central BH is suiTounded by an accretion disk whose radia- 
tion, with luminosity Lj, photoionizes the BLR, modelled as a 
spherical shell located at distance Rblr from the BH. Following 
GhiselUni & Tavecchio (2009), we set Rblr = ^^"l^als 
This relation provides a good approximation of the most recent 
results of the reverberation mapping studies (e.g. Kaspi et al. 
2007, Bentz et al. 2009). We suppose that the BLR clouds inter- 
cept and reprocess (mainly into emission lines) a fraction ^blr 
of Ld- As discussed in Tavecchio & Ghisellini (2008) a rather 
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good approximation for the BLR radiation field as seen in the 
comoving frame is a black body peaked at Vg^^j^ 3 x lO^^F Hz. 

Since the UVOT data-points probably trace the direct disk 
emission we fix by reproducing the UVOT fluxes with a 
black body. Assuming that the peak is in correspondence with 
the UVOT filter at the highest frequency (UVW2), a lower limit 
for the disk luminosity is Ld = 5 x 10"*^ erg s"' (Fig. 3, black 
short dashed line). This luminosity is exactly ten times larger 
than that estimated by Tanaka et al. (2011) by using the BLR 
total luminosity (and assuming ^blr ~ 0.1), in turn estimated by 
the luminosity of the H/3 line in the SDSS spectrum (Fan et al. 
2006). Considering the difference between the flux at the epoch 
of the SDSS (Jan 2008) and UVOT measures discussed in §2 
the discrepancy is reduced by a factor of 2. A difference by 
a factor of 5 between the two estimates could be explained by 
assuming ^blr ~ 0.02. Setting = 5 X 10"*^ erg s"', we have 
Rblr = 7 x 10'^ cm. 

Outside the BLR a dusty torus intercepts and re-emits part 
of the central disk emission. The radiation field of the torus is 
modeled as a black body with temperature Tjr = 1.2 x 10^* K 
with total luminosity Lir = 10'**' erg s"' (black long dashed line 
on Fig. 3) filling a volume that, for simplicity, is approximated 
as a spherical shell with radius Rm = 7 x 10"^ cm. 

In all the cases we model the compact emission region as a 
sphere with radius (the subscript "b" marks all the physical 
quantities related to the blob), moving with bulk Lorentz factor 
Fb, filled with uniform and tangled magnetic field B^. We assume 
that relativistic electrons follow a smoothed broken power law 
energy distribution with normalization K and slope ni and n2 
below and above the break at the energy ypineC^. 

We assume a conical jet propagating from the BH vicinity, 
with semi-aperture angle = 0.1 rad. A spherical region at 
distance d and radius R - (fid moving along the jet with bulk 
Lorentz factor F, carries tangled magnetic field (with uniform 
intensity B) and relativistic electrons with a distribution of the 
same functional form as in the case of the blob. 

In case (A) we consider only the emission from the compact 
blob. In case (B) and (C) we consider the emission from both 
regions. 

Both regions will emit synchrotron and IC radiation. For the 
blob and the jet components we consider as seed photons for the 
IC scattering the locally-produced synchrotron photons (SSC), 
the thermal emission from the dusty torus (EC/IR) and, in case 

(B) , the synchrotron photons produced in the other region of the 
jet (we will call them the EC/J and EC/b component). In case 

(C) for the large region we consider also the photons reprocessed 
from the BLR (EC/BLR). The calculations have been performed 
adapting the code described in detail in Maraschi & Tavecchio 
(2003). 

Since by construction the contribution of the large region is 
always negligible above 20-30 GeV, for simplicity we neglect 
absorption for this component when considering the region in- 
side the BLR. For the blob we expect a very weak absorption 
due to the IR radiation field of the torus, starting to be important 
at energies higher than a few TeV We also neglect this effect in 
the spectra shown in Figs. 3 and 4. 
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Fig. 2. Sketch of the geometrical arrangement assumed in the 
model (not to scale). We consider a spherical BLR with radius 
7?BLR and a dusty torus at R\^. In all the cases we consider the 
emission from a small compact "blob" of radius R^, moving with 
Lorentz factor Fb. While in case A we suppose that the blob is 
responsible for the entire SED, in case B and C we also consider 
the emission from a "standard" spherical emission with radius 
R equal to the cross sectional size of a conical jet with semi- 
aperture angle (p, moving with bulk Lorentz factor F located out- 
side (B) or inside (C) the BLR. Each region is characterized by 
different values of the physical parameters, such as the magnetic 
field, electron density and energies. See text for more details. 



3.3. Results 

We first consider (case A) the possibility that the very compact 
blob emitting TeV photons is responsible for the entire SED, 
from IR up to TeV energies (see Fig. 3). As a consequence of 
the large compactness implied by the short variability timescale, 
the system naturally produces a powerful SSC component peak- 



ing in the X-ray band. As detailed in the Appendix A, the short 
variability timescale and the condition that the SSC compo- 
nent lies below the observed X-ray spectrum constrains the ratio 
5b/<Jb ~ lO"*^ G. Therefore either we adopt a "standard" value 
of the magnetic field {B\, = 0.1 - 1 G) using extremely large 
Doppler factors, 6b > 50 or, conversely, we fix Sb to smaller 
values 5b ~ 10 with extremely low magnetic fields, B^ ~ 10""* 



4 



F. Tavecchio et al.: On the origin of gamma-ray emission of PKS 1222+216 



— 8 I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




Log u [Hz] 

Fig. 3. As in Fig. 1 but with the model considering only a 
compact emission region located beyond the BLR (see text for 
details). Black short and long dashed curve show the assumed 
emission from the accretion disk and dusty torus. For clarity the 
XRT spectrum of May 29 is not reported. 



G. In view of the possibility that the blob is the result of recon- 
nection events in the jet, expected to produce rather fast outflows 
(Giannios et al. 2009), in our model we adopt the first possibility, 
using 6b -75. 

As commonly found in FSRQs (e.g. the discussion in Celotti 
& Ghisellini 2008), the condition that the the low energy tail 
of the EC component reproduces the X-ray spectram implies 
that the value of -ytnin is constrained to be close to unity. As a 
consequence the large number of particles at low energy pushes 
the total power demand of the region (calculated assuming one 
cold proton per electron) to extremely large values, close to 
f jet = 3 X 10"^^ erg s"' (see Table[T]i. It would be possible to de- 
crease the requested power allowing larger values of the Doppler 
factor ^b- However, even with an extreme value of 6b = 100 the 
power remains quite large. Such large jet powers are not uncom- 
mon among FSRQs in bright states (e.g. Ghisellini et al. 2010; 
see also Bonnoli et al. 2011 for the specific case of 3C454.3, 
which has a similar L^) and therefore we cannot exclude this 
possibility for PKS 1222+216. 

In Fig. 4 we present the results for the two-zone model in the 
two cases discussed above. In case B (left) we assume that both 
regions are cospatial and located outside the BLR, at distances 
d > /?blr; in case C (right) the large region is inside the BLR, 
while the blob is kept at the same distance as before. For the 
purposes of the calculation the exact distance of the two regions 
is not important, since the radiation field of the BLR and the 
torus are uniform within the corresponding radii. 

As demonstrated in Appendix B, in case (B) we find that 
the radiative coupling between the two regions is less impor- 
tant than the EC flux using IR photons as seeds for the IC emis- 
sion. Therefore the EC/J and EC/b components can be safely 
neglected. 

Since, in both cases, we are dealing with two separate re- 
gions, the number of parameters is relatively large and therefore 
we have some freedom when selecting their values. However, 



the choice is not arbitrary: in the emission models shown in Fig. 
4 we were guided by the following criteria. 

Compact blob — We assume that it dominates the y-ray 
emission but it provides a small contribution in the X-ray band 
(but see below). This allows to fix yp, to put a lower limit to 
Tmax, and to impose a lower limit on the minimum Lorentz fac- 
tor of the emitting electrons (thus decreasing the total power re- 
quest). Similarly to case A (see discussion above) we assume a 
large bulk Lorentz factor to keep the SSC emission below the ob- 
served X-ray luminosity. Finally the magnetic field and the elec- 
tron density are determined by the level of the y-ray emission 
and by the condition that the synchrotron emission lies below 
the observed optical flux level. 

Large region — For the IC component we are guided by the 
X-ray spectrum and by the observed flux. Moreover, we assume 
that its MeV-GeV emission contributes substantially to the total 
observed flux, although we assume that it lies a factor x 2 below 
the blob emission since the the LAT and the MAGIC spectra 
suggest that they describe a unique component. The synchrotron 
peak is constrained by the new IR data. For case C (emission 
inside the BLR) we tried to use values of the physical parameters 
commonly inferred in FSRQs (e.g. Ghisellini et al. 2009). 

Concluding, while parameters cannot be uniquely deter- 
mined, their values are reasonably constrained by the conditions 
posed by the SED and the observed variability. The parameters 
used for the SEDs shown in Fig. 3-4 are reported in Table [1] 
Doppler factors are calculated assuming that the jet is observed 
at an angle 0y = 0.7° from the axis. This relatively small an- 
gle is required in order to allow the assumed large value of the 
Doppler factor of the blob, 6b = 75. The last three columns re- 
port the power caiiied by cold protons, relativistic electrons and 
magnetic field. For the proton component we assume one pro- 
ton per relativistic electron. For cases B and C the power carried 
by the jet is in both cases exceedingly larger than that of the 
blob. This is somewhat different than the case of similar models 
studied for the BL Lac PKS 2155-304 in GhiselHni & Tavecchio 
(2008), where the power inferred for the blob was larger by more 
than one order of magnitude than that of the jet component. One 
reason for the difference is that the large external energy density 
makes the emission from the blob much more efficient than in the 
case of PKS 2155-304 (for which the dominant radiation field 
was that of the jet), thus reducing the required number of par- 
ticles. Moreover, the total number of protons is also minimized 
by the fairly large minimum Lorentz factor of the emitting elec- 
trons. We also note that in all cases, both for the blob and the 
jet, the resulting intensity of the magnetic field is rather low. In 
turn this also implies that the associated Poynting flux is neg- 
ligible compared to the power caiiied by the jet in the form of 
electrons and protons. The jet is therefore strongly matter domi- 
nated. This conclusion is robust, being the result of the the very 
large observed ratio between the high energy (IC) and low en- 
ergy (synchrotron) bumps. 

4. Discussion and Conclusions 

The detection of intense, hard and rapidly variable VHE y-ray 
emission from FSRQs, with thermal emission from the accretion 
disk and broad emission lines, poses a puzzle. The observed very 
rapid variability constraints the size of the emitting region to be 
small, suggesting that its location in the jet is inside the BLR. But 
in this case the high energy photons interact with the emission 
line photons targets, and get absorbed. 

The solution of this puzzle necessarily requires the existence 
of a compact region located outside the BLR, whose size is 
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Fig. 4. As in Fig. [3] but with the model corresponding to emission from the compact region and the standard jet. The blue short 
dashed line shows the emission from the compact region, while the red long dashed line reports the emission from the large region 
of the jet. The solid black line is the sum of the two. The left panel shows the case in which both emission regions are located outside 
the BLR (case B); the right panel represents the case in which the large region is within the BLR (case C). 



smaller than the expected cross sectional radius of the jet at that 
distance. This scenario follows similar ideas proposed to explain 
the ultra-fast VHE emission observed in the BL Lac objects PKS 
2155-304 and Mkn 501. However, the case of FSRQs is more 
constraining, since, to avoid the strong absorption by the BLR 
photons, the blob must be located at large distances, > 0.1 pc. 
In BL Lacs, instead, the relatively "clean" environment is rather 
transparent to y-ray photons (at least up to the TeV band, above 
which the thermal radiation from a dusty torus could affect the 
emission, e.g. Celotti et al. 1998) and the emission region could 
reside much closer to the central engine. 

As in the case of BL Lacs, also for FSRQs such compact 
blobs could originate from reconnection events which could nat- 
urally produce compact regions of fastly moving plasma inside 
the jet (the so called "jet in a jet" scenario, Giannios et al. 2009, 
2010, Nalewajko et al. 201 1). However, with respect to the case 
of BL Lac objects, in which the emission and the electron cool- 
ing are dominated by the SSC process, the denser external ra- 
diation field is expected to dominate the radiation energy den- 
sity, leading to stronger radiative losses, lower average particle 
energies and, consequently, an emission peaking at GeV ener- 
gies. The derived total power carried by the "blob" is much 
less than the total jet power, supporting the energetic sustain- 
ability of the process, and the possibility that there are, at any 
given time, several active blobs, pointing in slightly different 
directions. This would increase the otherwise small probabil- 
ity to have very small viewing angles (see also discussions in 
Ghisellini et al. 2009; Giannios et al. 2010). It is conceivable 
that the events leading to the production of these fast blobs can 
occur also at smaller distances from the central BH, even within 
the BLR. In this case the dense radiation field would prevent the 
escape of the VHE photons and the entire power emitted above 
few tens of GeV would be eventually reprocessed at lower en- 
ergies. A non-negligible magnetic field within the BLR would 
spread the directions of the emitting pairs into a broad cone, thus 



decreasing the observed flux and the chances of detection of this 
reprocessed radiation. Although this scenario is attractive, it is 
based on the assumption that the plasma in the emission region 
is magnetically dominated, with magnetization parameter (i.e. 
Poynting to kinetic flux ratio) of the order of 10-100. On the 
contrary, our models show that the flow at the distances where 
radiation is produced is characterized by a rather low value of the 
magnetic flux compared to that carried by electrons or protons. 

The severe radiative losses of the electrons in the disk and 
BLR radiation fields seem to rule out also the alternative pos- 
sibility that the rapidly varying VHE emission originates in ex- 
tremely collimated beams of ultra-relativistic electron, as envis- 
aged in the "needle" scenario of Ghisellini et al. (2009). In fact, 
due to the strong dependence of the maximum electron Lorentz 
factor resulting from the competition between the magneto- 
centrifugal acceleration in the black hole magnetosphere (e.g., 
Rieger & Mannheim 2000) and the strong IC cooling on the ra- 
diation energy density, ymax (Osmanov et al. 2007), the 
presence of a luminous disk and BLR will result in rather low en- 
ergy of the electrons, emitting at energies below the GeV band. 

Another possibility to explain the existence of a very 
small region at large distance is by the focusing of the jet 
through strong recollimation (e.g. Komissarov & Falle 1997, 
Sokolov et al. 2004 Stawarz et al. 2006, Nalewajko & Sikora 
2008, Bromberg & Levinson 2009). In particular, Bromberg & 
Levinson (2009) showed that if the shock resulting from the rec- 
ollimation of the jet by the interaction with an external medium 
is efficient in converting the dissipated energy into radiation, 
the reduced post-shock pressure allows the jet to be strongly 
"squeezed", reaching a minimum size close to the point where 
the recollimation shock crosses the jet axis. At this point reflec- 
tion (internal) shocks are expected to form (e.g., Komissarov 
& Falle 1997), and the shortest variability time should be con- 
strained by the size of these structures, of the order of the ra- 
dius of the "nozzle" of the jet. Under the assumption that the 
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75 


50 
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8.2x10"" 


Jet (in) 


1 


300 


6x 10* 


2 


3.45 


0.75 


1.8 X 10" 


1.5 X lO"" 


19.7 


10 


81 


0.34 


6.3 x 10"- 


Blob 
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700 


4x 10^ 


2.2 


3.4 


0.1 


6x 10^ 


7.2 X 10'" 


75 


50 


10.6 


0.4 


6.5x10-' 


Jet (out) 


3 


2 X 10' 


6x 10" 


2 


4.1 


0.09 


3x 10" 


1.1 X 10" 


19.7 


10 


17.7 


0.3 


5 X 10-- 


Blob 


100 


900 


4 X 10' 


2.2 


3.6 


0.18 


10^ 


6.2 X 10'" 


75 


50 


3.9 


0.8 


1.5x10-" 



Table 1. Input parameters for the emission models shown in Fig. 3 (first row) and Fig. 4 (other rows); "in" and "out" means inside 
or outside the BLR. The parameters are: the minimum (ymin), break (yp) and maximum (ymax) Lorentz factor and the low energy 
(n\) and the high energy (n2) slope of the smoothed power law electron energy distribution, the magnetic field B, the normalization 
of the electron distribution, K, the radius of the emission region, R, the Doppler factor 6 and the corresponding bulk Lorentz factor 
F. Doppler factors are calculated assuming that the observer lies at an angle dy = 0.7° from the jet axis. The last three columns 
report the power canied by (cold) protons (assuming one proton per emitting electron), relativistic electrons and magnetic field in 
units of 10'*^ erg s-'. 



external confining medium is associated with matter in the BLR, 
Bromberg & Levinson (2009) estimate that the re-collimation 
shock reaches the jet axis at a distance d* of the order of d* ^ 
2.5 Lj 46 (Rblr/0.1 pc)-^ pc where Lj is the jet power. Their sim- 
ulations show that at this point if (/) the cooling operates at a rate 
larger than that of the dynamical timescale and (//) the radiative 
conversion efficiency is sufficiently large (30 percents of the bulk 
luminosity converted into radiation), the cross sectional radius 
of the jet in the focusing point can be as small as R ~ lO-^^ii*. 
Condition (/) is easily fulfilled in FSRQs, even outside the BLR. 
Nalewajko & Sikora (2009) showed that the efficiency of inter- 
nal energy production in recollimation shocks can be larger than 
30% provided the product of the jet bulk Lorentz factor and the 
aperture angle is larger than ^ 3. Therefore, provided that the 
further conversion of this dissipated energy into radiation can 
reach large efficiencies, also condition (//) could be satisfied. 
Considering the rough consistency of these numbers with the 
results of our radiative models and the difficulties faced by the 
"jet in a jet" or "needles" models discussed above, we believe 
that such a scenario is rather promising and worth to be better 
investigated. 

Acknowledgements. We are grateful to an anonymous referee for rather con- 
structive comments that leaded to substantially improve the paper We thank 
Amir Levinson for discussions. 

References 

Abdo, A. A., et al. 2010, ApJ, 722, 520 

Abdo, A. A., et al. 2011, Science, 331, 739 

Aharonian, R, et al. 2007, ApJ, 664, L71 

Albert, J., et al. 2007, ApJ, 669, 862 

Albert J. et al., 2008, Science 320, 1752 

Aleksic J., et al., 201 la, A&A, 530, 4 

Aleksic J., et al., 201 lb, ApJ, 730, L8 

Aloy M. A., Janka, H.-T., Mueller, E. 2005, A&A, 436, 273 

Bednai-ek, W., & Idee, W. 201 1, MNRAS, 551 

Bentz M. C, et al., 2009, ApJ, 705, 199 

Blazejowski M., Sikora M., Moderski R., Madejski G. M., 2000, ApJ, 545, 107 
Bloom S. D., Marscher A. R, 1996, ApJ, 461, 657 

Bonnoh G., Ghiselhni G., Foschini L., Tavecchio R, Ghirlanda G., 2011, 

MNRAS, 410, 368 
Bromberg O., Levinson A., 2009, ApJ, 699, 1274 
Celotti A., Fabian A. C, Rees M. J., 1998, MNRAS, 293, 239 
Celotti, A., & Ghisellini, G. 2008, MNRAS, 385, 283 
Chonis T. S., Gaskell C. M., 2008, AJ, 135, 264 
de Angehs A., Mansutti O., Persic M., 2008, Nuovo Cimento, 31, 187 
Decarh R., Dotti M., Treves A., 201 1, MNRAS, 413, 39 
Dermer C. D., Finke J. D., Krug H., Bottcher M., 2009, ApJ, 692, 32 
Dommguez A., et al., 2011, MNRAS, 410, 2556 
Fan Z., Cao X., Gu M., 2006, ApJ, 646, 8 

Foschini L., Ghisellini G., Tavecchio R, BonnoH G,, Stamen'a A., 2011, A&A, 
530, 77 



Ghisellini, G., Celotti, A., Fossati, G., Maraschi, L., & Comastri, A. 1998, 

MNRAS, 301,451 
Ghisellini G., Tavecchio R, 2008, MNRAS, 386, L28 

Ghisellini, G., Tavecchio, R, Bodo, G., & Celotti, A. 2009, MNRAS, 393, L16 
Ghisellini G., Tavecchio R, 2009, MNRAS, 397, 985 
Ghisellini G., et al., 2011, MNRAS, 411, 901 

Ghisellini G., Roschini L., Volonteri M., Ghirlanda G., Haardt R, Burlon D., 

Tavecchio R, 2009, MNRAS, 399, L24 
Giannios D., Uzdensky D. A., Begelman M. C, 2009, MNRAS, 395, L29 
Giannios D., Uzdensky D. A., Begelman M. C, 2010, MNRAS, 402, 1649 
Liu H. T., Bai J. M., 2006, ApJ, 653, 1089 
Hinton J. A., Hofmann W., 2009, ARA&A, 47, 523 
Jarvis M. J., McLure R. J., 2006, MNRAS, 369, 182 

Kaspi S., Brandt W. N., Maoz D., Netzer H., Schneider D. R, Shemmer O., 2007, 
ApJ, 659, 997 

Komissarov, S. S., & Palle, S. A. E. G. 1997, MNRAS, 288, 833 
Komissarov, S. S., & Lyutikov, M. 2011, MNRAS, 518 

Kubo, H., Takahashi, T, Madejski, G., Tashiro, M., Makino, R, Inoue, S., & 

Takahara, R 1998, ApJ, 504, 693 
Malmrose M., Marscher A., Jorstad S., Nikutta R., Elitzur M., 2011, ApJ, 732, 

116 

Maraschi L., Ghisellini G., Celotti A., 1992, ApJ, 397, L5 
Maraschi L., Tavecchio R, 2003, ApJ, 593, 667 
Marscher A. R, et al., 2008, Nature, 452, 966 

Marscher, A. P., & Jorstad, S. G. 2010, in Proc. Rernii Meets Jansky 

larXiv:1005.5551 1 
Nalewajko K., Sikora M., 2009, MNRAS, 392, 1205 

Nalewajko, K., Giannios, D., Begelman, M. C, Uzdensky, D. A., & Sikora, M. 

2011, MNRAS, 413, 333 
Osmanov, Z., Rogava, A., & Bodo, G. 2007, A&A, 470, 395 
Osterbrock, D. E., & Pogge, R. W. 1987, ApJ, 323, 108 
Padovani P, Giommi P, 1995, ApJ, 444, 567 
Poutanen J., Stem B., 2010, ApJ, 717, LI 18 
Rieger, R M., & Mannheim, K. 2000, A&A, 353, 473 

Shields, G. A. 1978, in Proc. of the Pittsburgh Conf on BL Lac objects, Ed. 

A.M. Wolfe, p. 275 
Sikora, M., Begelman, M. C, & Rees, M. J. 1994, ApJ, 421, 153 
Sikora M., Moderski R., Madejski G. M., 2008, ApJ, 675, 71 
Sokolov, A., Marscher, A. P, & McHardy, I. M. 2004, ApJ, 613, 725 
Spada, M., Ghisellini, G., Lazzati, D., & Celotti, A. 2001, MNRAS, 325, 1559 
Stawarz L., Aharonian F., Kataoka J., Ostrowski M., Siemiginowska A., Sikora 

M., 2006, MNRAS, 370, 981 
Stem, B. E., & Poutanen, J. 201 1, MNRAS, submitted (arXiv:l 105.2762t 
Tagliaferri, G., et al. 2008, ApJ, 679, 1029 
Tanaka Y. T, et al., 2011, ApJ, 733, 19 
Tavani, M., et al. 2011, Science, 331, 736 
Tavecchio R, Maraschi L., Ghiselhni G., 1998, ApJ, 509, 608 
Tavecchio, R, et al. 2000, ApJ, 543, 535 
Tavecchio, R, et al. 2002, ApJ, 575, 137 
Tavecchio R, Ghiselhni G., 2008, MNRAS, 386, 945 
Tavecchio R, Mazin D., 2009, MNRAS, 392, L40 

Tavecchio, R, Ghisellini, G., Bonnoh, G., & Ghirlanda, G. 2010, MNRAS, 405, 
L94 

Tomikoski M., et al., 1996, A&AS, 116, 157 
Vittorini, V., et al. 2011, ApJ, 732, L22 

Wagner, S., & Behera, B. 2010, in 10th HEAD Meeting, Hawaii (BAAS, 42, 2, 
07.05) 



7 



F. Tavecchio et al.: On the origin of gamma-ray emission of PKS 1222+216 



Appendix A: Limit to B and 6 from the constraints 
on the SSC component 

One can derive an approximate relation between the magnetic 
field and the Doppler factor (assuming for simphcity 6^, = Ft) 
starting from the approximate expressions for the EC and SSC 
luminosities: 



and: 



(A.1) 



where now primed quantities are measured in the jet frame. As 
above, this expression can be rewritten as: 



^EC/IR _ ilR / ^1 
^EC/b \^Ir/ \^iel 



5b \^ 1 



(B.5) 



Again, since Ljs. is comparable to or larger than the observed Lb , 
Rm~lOx R], Tb = 50, Fj = 10, 6] ~ 20, LEc/m/LEc/b » 1- 

Therefore both EC/J and EC/b components can be safely ne- 
glected. 



Lec = TO-TcU^nTl <r'>nV6t 



(A.2) 



where U'^^^^^ and Uext^l are the energy density of the syn- 
chrotron and external radiation in the blob frame, V = 4/3;r/?^ 
is the blob volume, n is the electron density and < > is the 
averaged value of the squared of the electron Lorentz factor. 

Using these equations imposing that Lec - Ly, Lssc /Lec < 
Lx/Ly and Rb - c/var^b. with some manipulations one can de- 
rive the following relation between the magnetic field and the 
Doppler factor (assuming for simplicity 6b = Fb) containing only 
observed quantities: 



Bb < 





LxV 


Ly \ 


Ly) 


r specific 


case 



1/2 



(A.3) 



I = 600 s, Ly ^ 3 X erg s-\ 
Lx/Ly < 10 and f/ext = Um = 5 X 10"* erg cm~^ we derive 
fib - 10-4(5b/ 10)5 G. 



Appendix B: Importance of the blob/jet radiative 
interplay 

The ratio of the EC/IR and EC/J (inverse Compton produced in 
the blob scattering photons of the jet) luminosities is given by the 
ratio of the energy densities of the corresponding seed photons 
in the blob rest frame: 

^EC/lR ,„ ^, 

1 ~ IF ^ - ^ 

i-EC/J 

where primed quantities are measured in the blob frame. The 
two energy densities can be expressed as: 

U'lR - U', = (B.2) 



where "b" and "J" refers to quantities related to the blob and the 
jet, respectively, and F^i = FbFj(l - pbPi) is the relative bulk 
Lorentz factor between the jet flow and the blob. Therefore we 
finaUy find: 

(B.3) 

Leg/ J Lj \RiRj \^xe\l 

Since Lir is comparable to the observed Lj (see Fig. 1), 
/?iR ~ 10 X Fb = 50, Fj = 10, 5j ~ 20, the ratio is largely in 
favor of the EC/IR component, Lec/ir/Lec/j ~ 10*. 

Analogously, one can evaluate the ratio of the EC/IR and 
EC/b (inverse Compton produced in the jet scattering photons of 
the blob) luminosities: 

■f-EC/b f^b 
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